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THE INVISIBLE MONSTER HAS TWO FACES: 
OBSERVATIONS OF e AURIGAE WITH THE HERSCHEL SPACE OBSERVATORY 
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ABSTRACT 

We present Herschel Space Observatory photometric observations of the unique, long-period eclipsing bi- 
nary star e Aurigae. Its extended spectral energy distribution is consistent with our previously published cool 
(550 K) dust disk model. We also present an archival infrared spectral energy distribution of the side of the 
disk facing the bright F-type star in the binary, which is consistent with a warmer (1 150 K) disk model. The 
lack of strong molecular emission features in the Herschel bands suggests that the disk has a low gas-to-dust 
ratio. The spectral energy distribution and Herschel images imply that the 250 GHz radio detection reported 
by Altenhoff et al. is likely contaminated by infrared-bright, extended background emission associated with a 
nearby nebular region and should be considered an upper limit to the true flux density of e Aur. 

Subject headings: stars: AGB and post-AGB — binaries: eclipsing — circumstellar matter — stars: individual 
(Epsilon Aurigae) 



1. INTRODUCTION 



In iHoard et"aTJ (120101) (henceforth, Paper I), we provide a 
summary of previous results on e Aurigae, the eclipsing bi- 
nary star with the longest known orbital period (27. 1 yr), and 
introduce an empirical model for this system. Our model 
was constrained by a spectral energy distribution (SED) span- 
ning three orders of magnitude in wavelength, from the 
far-ultraviolet to the mid-infrared, assembled from new and 
archival ground- and space-based observations. We estab- 
lished the nature of e Aur as a large (R sw 135 Rq), but low 
mass (M ss 2 Mq), post-asymptotic giant branch F-type star 
orbiting a B5 main sequence star surrounded by a large (R « 4 
AU), cool (550 K) dust disk. Subsequent to Paper I, a num- 
ber of authors have further illuminated the e Aur disk (e.g., 
Stencel et al]|2011l lBudaill2011l IChadimaet aTJl2011l iTakeutil 
20111) . highlighting common threads pointing to agreement 
with the low mass scenario and exi stence of discrete optically 
thick and thin disk components. Kloppenborg e t al.l (1201 Oh 
even obtained near-infrared interferometric images of the disk 
passing across the face of the F star during eclipse ingress, 
which have enabled better constraints on the characteristics 
of the optically thick disk component. 

The Paper I SED was constructed primarily from data ob- 
tained close in time to the latest (2009-2011) or previous 
(1982-1984) primary eclipse of e Aur. As such, it repre- 
sents a view of the unilluminated "back" side of the disk rel- 
ative to the F star. We now present an expansion of the Pa- 
per I SED through another order of magnitude of wavelength, 
via new observation s with the Herschel Space Observatorj^ 
(Pilbra tt et alJl20IO). We also explore the properties of the 
"front" side of the disk (i.e., the side facing - and directly il- 
luminated by - the F star) via archival observations obtained 
at orbital phases around disk opposition. 



1 Spitzer Science Center, California Institute of Technology, Pasadena, 
CA 91125 

2 Department of Physics and Astronomy, University of Denver, Denver, 
CO 80208 

3 NASA Ames Research Center, Moffett Field, CA 94035 

4 Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and important participation 
from NASA. 



2. OBSERVATIONS 

We observed e Aur with the Photodetector Arra 



Camera and Spectrometer (PACS; Pogl itsch et alj 1201 
and Spectral and Photometric Imaging Receiver (SPIRE; 
iGriffin et alJ2010|) on Herschel (Open Time Cycle- 1 program 
"OTl_hoard_l"). These observations were accomplished on 
9-10 Sep 2011 UT (PACS) and 21 Sep 2011 UT (SPIRE), 
shortly after the end of the optical eclipse in early July 201 1 
dStencel et alj|201 lh . The observations consist of 2, 2, and 
4 "mini-scan maps" of 10 repetitions in each of the PACS B 
(Bl; 70 fxm), G (B2; 1 10 /im), and R (160 /im) bands, respec- 
tively, with each pair at cross-scan angles of 70 and 110 de- 
grees, and 8 "small map" repetitions at 250, 350, and 500 /im 
for SPIRE. 

We first fit each image using a Gaussian profile to precisely 
determine the source position. We then integrated around this 
position within apertures of 10, 15, 20, and 22 arcsec for the 
70, 110, 160, and 250 /im images, respectively (e Aur is not 
detected in the 350 and 500 /im images - see below). The 
aperture photometry was performed on Level 2 data products 
(reduction pipeline v7.2.0 for PACS and v7.3.0 for SPIRE) us- 
ing a custom photometry routine written by D. Ladjal that is 
based on the DAOPHOT photometry package but with error 
calculations adapted to bolometers. We applied corrections 
for both aperture (from the online PACS and SPIRE photom- 
etry cookbooks) and color (based on the slope of the photo- 
metric points as a function of wavelength). 

Figure Q] shows images from each Herschel band, and the 
photometry is listed in Table [T] The photometric uncertain- 
ties are quadrature sums of the absolute measurement errors 
(0.09%, 0.3%, 1.8%, 21% for 70, 110, 160, 250 /im, respec- 
tively) and the instrumental and calibration errors (5% for 
PACS-70, -110; 10% for PACS-160; 15% for SPIRE). As a 
check of the Herschel photometry, we note that the PACS- 
70 value differs by only 1% from the Spitzer MIPS-70 value 
from Paper I (listed here in Table [TJ, well within their mutual 
la uncertainties (5% and 15%, respectively). 

The upper limits in Table[T]at 350 and 500 /im are 3cr r.m.s. 
values for each map (in mJy/beam). Inspection of the images 
suggests that these non-detections can be blamed, at least in 
part, on the presence of infrared-bright structured background 
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FIG. 1. — Herschel Level 2 pipeline images of e Aur consisting of combined PACS cross-scan images at 70, 110, and 160 /im (top row) and SPIRE maps at 
250, 350, and 500 fim (bottom row; wavelengths increasing from left to right in both rows). The SPIRE images are shown with a logarithmic grayscale stretch, 
while the PACS images are linear. The SPIRE images are 6.4 arcmin wide by 5.9 arcmin tall (the circles marking the location of e Aur have a diameter of 40 
arcsec); the PACS images have the same angular scale as the SPIRE images. All of the images are oriented with north up and east to the left, rectified to the world 
coordinate solution of the PACS 70 fim image. 

emission, rather than simply faintness of e Aur at these wave- 
lengths. This structured background is likely related to the 
know n dark cloud nebular region TGU HI 105 (Dobashi et al. 
2005) centered -11 arcmin northeast of e Aur and extending 
southward - this is visible as the very bright extended emis- 
sion along the east edge of the SPIRE maps. The point source 
visible in the PACS images w 90 arcsec southeast of e Aur 
(J2000 coordinates of 05:02:03.6, +43:48: 15.9 measured from 
the PACS-160 image) is likely an embedded star or proto- 
star, since it is not visible in optical (DSS) or near-infrared 
(2MASS) imaged 

3. DISCUSSION 

3.1. Orbital Geometry 

The e Aur binar y orbit has a large eccentricity (0.23; 
IStefanik eT al. 2010). Consequently, times of primary eclipse 
and opposition phase of the disk do not occur exactly 
0.5P nr h apart, as the y would for a circular orbit. Instead, 
Stefan ik et alJ ((2010) showed that the red-to-blue crossing of 
the F star radial velocity curve (which is coincident with the 
photometric eclipse mid-point) occurs « 0.1P or b after the pe- 
riastron passage of the disk around the F star, and the blue-to- 
red crossin g occurs another » 0.6P or b later (e.g., see Figures 
1 and 4 in Stef anik et al.l l2010). Thus, when using a linear or- 
bital ephemeris with zeropoint at mid-eclipse, the disk oppo- 
sition occurs at <fi s=s 0.6. To determine when the data were ob- 
tained relative to the viewing geometry of the binary, we use 
a linear ephemeris constructed by combining the best orbital 

5 It is also not visible in the Wide-Field Infrared Survey Explorer (WISE) 
Preliminary Release Atlas images at 3.4 and 4.6 fim. A faint source is possi- 
bly present at this position in the 12 and 22 fim WISE images, but is not listed 
in the photometry catalog, indicating that, if real, it has S/N<7 in both bands 
and is fainter than ss 1 and fa 6 mjy, respectively. 



period found by IStefanik et al.l (120101) with the latest eclipse 
midpoint date found bv lStencel et al.1 (1201 ll) (see Table [TJ. 

3.2. Cool Side of the Disk 

Figure [2] shows an SED of e Aur derived from the one pre- 
sented in Paper I, but starting in the near-infrared at 1 /im. We 
increased the long wavelength limit compared to Paper I by 
including both the new Hers chel photometry and th e 250 GHz 
radio detection reported by lAltenhoff et al.l dl994l) . We also 
differentiated the data listed in Table 1 of Paper I into two 
groups that were measured at <fi w 0.85-0. 15 (as in Paper I, the 
SED ex cludes data obtained during eclipse, at <f> =0.97-0.03; 
e.g., see IStencel et al.ll20Tll) and at <j) wO.5-0.7. The former 
data correspond to the cool back side of the disk (facing away 
from the F star) during pre- and post-eclipse, while the latter 
were obtained around opposition phase and correspond to the 
warm front side of the dis k (directly facing the F star). The 
back side data exclude the lTaranova & S henavrin (200JJ) and 
Midcourse Space Experiment (MSX) photometry from Paper I 
(which were obtained at <p ~ 0.6). The exceptions to this divi- 
sion of the data are the near-infrared photometric points from 
Paper I, all of which were obtained at orbital phases close to 
0.6. However, the disk contribution is less than a few hun- 
dredths of a percent of the total model shortward of K (less 
than 3% even in the warm - hence, bright - disk case; see be- 
low); thus, the near-infrared photometry is effectively phase- 
independent. 

Figure [2] also shows the model from Paper I, extended to 
1500 /im. This model utilizes a parametric description of a 
homogeneous, partially transparent dust disk that is charac- 
terized by a "transmissivity" (the degree of transparency of 
the disk) and an "emissivity" (the degree to which light radi- 
ated by dust grains inside the disk is visible from outside the 
disk, not to be confused with the standard emissivity e). Both 
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TABLE 1 
Photometry of e Aur 



Wavelength 
(pm) 


Instrument/ 
Band 


Observatory 


Date of Observation 
(JD-2400000) 


Orbital 
Phase" 


Flux Density 
(mJy) 


Reference 


71.44 


MIPS-70 


Spitzer 


53639, 53790 


Ri0.83 


497±74 


Paper I 


70 


PACS 


Herschel 


55814 


0.043 


492±25 


this work 


110 


PACS 


Herschel 


55814 


0.043 


254±13 


this work 


160 


PACS 


Herschel 


55814 


0.043 


105±25 


this work 


250 


SPIRE 


Herschel 


55826 


0.044 


57±15 


this work 


350 


SPIRE 


Herschel 


55826 


0.044 


<6() (3<r) 


this work 


500 


SPIRE 


Herschel 


55826 


0.044 


<57 (3o-) 


this work 


1200 


radio 


IRAM 30-m 


47710-47863 


raO.23 


9±2 


Altenhoff et al. (1994) 



a Calc ulated using the ep hemeris JD b s = ID 2 455390( 10) + 9896(1. 6)E (mid-eclipse date from lStencel et alj 
l2TjTil orbital period from Stefanik et al. 20lTj|). 
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FIG. 2. — Long wavelength SED of e Aur, showing measurements obtained at orbital phases correspondi ng to the coo l (pr e- and post-eclipse) and warm 
(opposition phase) sides of the dust dis k. The data are as foll ows: I, H, K/Ks from 2MASS (Paper I), Backman et al. 1 1984) and Taranova & Shenavrin 120011) . 
as well as the 250 GHz radio detection I Altenhoff et al. 1994) (unfilled squares). The data corresponding to the cool disk case (unfilled circles) include: Spitzer 
IRAC 3.5, 4.5 pm and MIPS 24, 70 pm (Paper I), ground-based L', M, N, Q flBackman et~aLl l 19841), and our new Herschel PACS and SPIRE (plotted as upper 
limits for SPIRE 350 and 500 (im). The Spitzer IRS and MIPS-SED spectroscopic data (Paper I) are plotted in grey. The data corresp onding to the warm disk 
case (filled circles) include: ground-based L, M (Taranova & Shenavrin 2001) and the MSX Bl, B2, A, C, D, E bands (Eean et al. 2003). The vertical bars mark 
the Herschel PACS and SPIRE bands. Two models are plotted: the cool case model (lower, red solid line) from Paper I is the sum of limb-darkened model FO 
(post-AGB) and B5 V spectra (upper and lower dotted lines, respectively), and a cool (550 K) blackbody disk (lower, red dashed line); the warm case model 
(upper, blue solid line) is the same but utilizes a 1 150 K disk (upper, blue dashed line). The inset panel is an expanded view of the 1-50 (im region, with the 
y-axis scaled by u (i.e., divided by the slope of the Rayleigh- Jeans tail of the blackbody function) to emphasize the difference between the cool and warm 
cases, which are essentially indistinguishable outside of this wavelength range. 

parameters are simple multiplicative scale factors applied to a 
single-temperature blackbody disk model. This model is ef- 
fectively a bulk average over a disk composed of discrete op- 
tically thick and thin regions, that minimizes the need to con- 
strain a larger number of parameters corresponding to multi- 
ple disk components. The Herschel photometry is completely 
consistent with the 550 K dust disk model from Paper I; how- 
ever, the radio point is 3.5rr brighter than the model. 

The e Aur dust disk almost cert ainly does co n tain a mix 
of opt i cally thick and thin regions; IStencel et all (1201 ll) and 
iBudail d2011h have discussed this, and Kloppenbo rg et alj 



(2010) have estimated likely geometric parameters for the op- 
tically thick component based on their interferometric imag- 
ing. Similarly, the actual dust disk should obey a radial tem- 
perature profile characterized by higher temperature closer to 
the central star, even if the observational signature of this is 
likely to be relatively minimal, as the disk is viewed at high 
inclination (i w 89°; Paper I), so its SED is dominated by 
the geometrically thick edge of the disk. As a first step to- 
ward a more "physical" disk model for e Aur (a full treatment 
is beyond the scope of this pa per), we modified t he detailed 
dust disk model described in Hoa rd et ail (|2009) and refer- 
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ences therein, which was developed for modeling dust disks 
around white dwarfs and cataclysmic variables, to calculate a 
two-component (optically thick and thin) model for the e Aur 
disk. 

Usin g the values derived in Paper I and Kloppenborg et al. 
(2010), we constrained the optically thick component to have 
outer radius of 3.8 AU, thickness of 0.76 AU, and outer edge 
temperature of 550 K, with an inclination of 89° at a distance 
of 625 pc. The inner edge was constrained to have a temper- 
ature of 1500 K (the assumed dust sublimation temperature). 
Given the 15,000 K B5 V central star described in Paper I, 
and assuming a radial temperature profile exponent of 0.56, 
this results in an inner disk radius of 0.55 AU, consistent with 
the discussion in Paper I concerning a possible "hole" in the 
dust disk where the central star sublimates the dust into gas. 
The slight increase in the radial temperature profile exponent 
compared to the expected value of 0.5 (corresponding to the 
surface of an opaque blackbody disk), suggests that the sur- 
face of the disk is not completely opaque down to some geo- 
metric depth at which the optical depth becomes large. This is 
comparable to the effect encapsulated in the "transmissivity" 
and "emissivity" parameters of the Paper I model. 

The optically thin disk component is modeled as a corona or 
atmosphere of hotter dust grains extending 0. 1 AU above and 
below the optically thick disk component. This component 
is assumed to be composed of 50-micron silicate (3 g cm" 3 ) 
grains, constrained to the same inner edge boundary condi- 
tions as the optically thick component, but having an outer 
edge temperature of 1000 K (corresponding to R = 1.9 AU). 
The total mass of dust in this component is 0.000 7 M ffi , 1% of 
the mass estimated by Kloppenbor g et all (120101) for the opti- 
cally thick disk component in e Aur. The combined contri- 
bution of these two model components produces a total SED 
almost indistinguishable from the Paper I cool disk model, 
differing by, at most, 2.8% (at 5 /im) and by less than 1% 
over 85% of the 1-1500 /im wavelength range. Although far 
from a comprehensive system model for e Aur (e.g., it lacks 
strong constraints on the additional degrees of freedom intro- 
duced by the increased physical complexity), the agreement 
between this calculation and the "bulk-averaged" model from 
Paper I provides some confidence that the latter is useful for 
making reasonable estimates of the disk parameters in e Aur. 

3.3. Warm Side of the Disk 

Using the warm side data plotted in Figure|2j we calculated 
a new version of the Paper I model to reproduce the SED of 
the front side of the e Aur disk (i.e., the side closer to, and 
directly illuminated by, the F star). This model is identical 
to the cool disk model shown here and in Paper I, except that 
the disk temperature is increased to 1150 K and the "emis- 
sivity" factor is set to 1.00 (down from 2.43 in the cool disk 
model). The disk temperature uncertainty, determined by not- 
ing the range over which the match between model and data 
beco mes noticeably worse, is ±50 K. The 250 GHz radio 
point ( Altenhoff et al. 1994) is still several a brighter than the 
warm disk model, which does not differ significantly from the 
cool disk model at very long wavelengths. We note in pass- 
ing that the relatively faint 12 and 20 /im MSX points in the 
warm case data - echoed by the low Q band point from the 
cool case data - might point to the presence of silicate ab- 
sorption features that were visible in the orbital cycle leading 
up to the 1982-1984 eclipse (when the photometric data were 
obtained) but not, apparently, visible in the IRS spectrum ob- 
tained prior to the latest eclipse. 



The higher temperature is an expected consequence of irra- 
diation from the F star. The lower "emissivity" might be un- 
derstood if increased radiation pressure from the F star on the 
front side of the disk compresses the disk material, effectively 
increasing its optical depth to the point at which the observ- 
able side of the disk at opposition is completely opaque. This 
opacity change could also account for the steepness of pho- 
tometric egr ess compared to in gress - see the light curve in 
Figure 2b of St encel et al.l (1201 ll) . For completeness, we note 
that an alternative scenario is that the "emissivity" could be 
> 1 if the warm disk emitting area is not completely opaque, 
but is smaller than the cool side emitting area. Exploring this 
possibility in a more quantitative fashion is precluded by the 
sparseness of the data at the longest wavelengths in the warm 
case SED. As described in Paper I, the disk "transmissivity", 
which is unchanged in the warm model compared to the cool 
model, is primarily important in the ultraviolet and has little 
effect on the model's match to the infrared data. 

4. DISCUSSION 

Takeutjl d201 ll) presents an analytic model for the disk in 
e Aur that accounts for both the geometry and rotation of the 
disk and irradiation by the F star. In Takeuti's model, the ini- 
tial system parameters are based on our empirical disk model 
from Paper I, but the disk temperature, rather than being de- 
fined, is a calculated quantity. For a particular assumed value 
for the heat capacity of the disk material that produces a tem- 
perature of «500 K on the back side of the disk (comparable 
to the empirical result from Paper I and this work), Takeuti's 
model predicts a disk front side temperature of «1200 K, in 
close agreement with our empirical result of 1 150 ± 50 K. 

In light of the Herschel extension to the e Aur SED, which 
closely follows the disk model from Paper I and the op- 
tically thick/thin model presented here, the Altenh off et alJ 
(1994) radio detection should be re-examined. Admittedly, 
we cannot rule out the possible existence of a very cold, or 
non-thermal, component that contributes significantly only at 
wavelengths longer than ~ 1000 /im. In addition, because the 
radio observation was made at an orbital phase close to disk 
quadrature, it is unclear whether it is better associated with 
the cool or warm side of the disk. The long wavelength end of 
the former is better constrained by data than that of the latter; 
thus, while the radio point seems incongruous in comparison 
to the rest of the cool disk SED, again, we cannot rule out 
that it originates from an unexpected component that is only 
visible on the warm side of the disk. However, we suggest 
that a more likely scenario is that the 250 GHz radio detec- 
tion (which had a half -power beam width of w 12 arcsec) is 
contaminated by - or potentially due solely to - the infrared- 
bright structured background in the vicinity of e Aur visible 
in t he Herschel SP I RE m aps (see Figure Q]). In that regard, 
the [Alt enhoff et al. ( 1994) value should be considered an ex- 
treme upper limit to the true flux density of e Aur at this long 
wavelength. Based on our model, we predict a much lower 
1200 /im flux density of « 2 mJy for e Aur. 

5. CONCLUSIONS 

A conclusion that can be drawn immediately from the ex- 
tended cool disk SED of e Aur in Figure [2] is that any poten- 
tially bright emission features from molecular species in the 
PACS and SPIRE bands are minimal, given the adherence to 
the 550 K disk model. This leads us to suggest that the gas- 
to-dust ratio in the e Aur disk is very low. The disk itself 
is likely of low mass even when the gas content is included 
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(<C 1 M Q , as suggested in Paper I and iKloppenborg et alj 
2010), and is perhaps analogo us to the transitional disk in 
HD 95881 dVerhoeffetalJl2010h . 

During the course of the several years between the end of 
primary eclipse (late 201 1), disk quadrature (circa 2019), and 
secondary eclipse (circa 2027), the side of the e Aur disk that 
is directly illuminated and heated by the F star is increasingly 
visible. Infrared and/or sub-mm observations during this in- 
terval, providing longitude by longitude characterization of 
the disk's azimuthal temperature profile, can uniquely address 
disk thermal gradients that are not feasibly observed at any 
other time until after the next eclipse. We urge observers to 
seize this rare opportunity. 

This work is based on observations made with Herschel, 



an ESA Cornerstone Mission with significant participation 
by NASA, and with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under a contract with NASA. Support 
for this work was provided by NASA through an award is- 
sued by JPL/Caltech. We used data products from the Two 
Micron All Sky Survey, a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis Cen- 
ter/Caltech, funded by NASA and NSF, and utilized the SIM- 
BAD database, operated at CDS, Strasbourg, France, and 
NASA's Astrophysics Data System. We thank Sean Carey for 
a helpful discussion about MSX photometry. 

Facilities: AAVSO, Herschel (PACS, SPIRE), IRAS, MSX, 
Spitzer (IRAC, IRS, MIPS) 
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